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ABSTRACT 

Advection-dominated, high-temperature, quasi-spherical accretion flow onto a com- 
pact object, recently considered by a number of authors, assume that the dissipation of 
turbulent energy of the flow heats the ions and that the dissipated energy is advected 
inward. It is suggested that the efficiency of conversion of accretion energy to radia- 
tion can be very much smaller than unity. However, it is likely that the flows have an 
equipartition magnetic field with the result that dissipation of magnetic energy at a rate 
comparable to that for the turbulence must occur by Ohmic heating. We argue that 
this heating occurs as a result of plasma instabilities and that the relevant instabilities 
are current driven in response to the strong electric fields parallel to the magnetic field. 
We argue further that these instabilities heat predominantly the electrons. We conclude 
that the efficiency of conversion of accretion energy to radiation can be much smaller 
than unity only for the unlikely condition that the Ohmic heating of the electrons is 
negligible. 

Subject headings: accretion, accretion disks — galaxies: active — plasmas — magnetic fields — 
stars: magnetic fields — X-rays: stars 



1. Introduction 



2. Accretion Flows with B Field 



Advection-dominated accretion flows have been in- 
tensely studied during the past several years (for ex- 
ample, Narayan and Yi 1994; Narayan and Yi 1995; 
Abramowicz ct al. 1995; Nakamura et al. 1996; 
Chakrabarti 1996). The basic dynamical equations 
for accretion disks including the advection of entropy 
were first disciissed by Paczyhski and Bisnovatyi- 
Kogan (1981) and Muchotrzeb and Paczyhski (1982). 
In contrast with the widely applied theory of thin ac- 
cretion disks of Shakura (1973) and Shakura and Sun- 
yaev (1973) where the disk material cools efficiently 
by local radiation of viscously generated energy, the 
advection-dominated accretion fiows of Narayan and 
Yi assume that the viscous dissipation heats the ions 
and that a constant fraction / of this dissipated en- 
ergy is advected inward and that the fraction 1 — / 
is locally radiated. The further assumption that the 
energy exchange between ions and electrons is by 
Coulomb scattering leads to conditions with the ion 
temperature Tj much larger than the electron temper- 
ature Te so that the cooling is inefficient. [A recent 
paper by Esin et al. (1996) treats advection domi- 
nated accretion flows assuming Tj = Tg.] The radia- 
tive efficiency, the power output in radiation divided 
by Mc^ (with M is the mass accretion rate), is found 
to be very small compared with unity. The advection- 
dominated accretion flows tend to be quasi-spherical 
and optically thin (except for cyclotron radiation as 
discussed below) with radial inflow speed Vr » —avK, 
azimuthal speed « const. ^ vk, ion thermal 
speed Cgi « const. ~ vk (Narayan and Yi 1995), 
where Vk = {GM/r)^ is the Kepler speed and a is the 
dimensionless viscosity parameter of Shakura (1973) 
usually assumed to be in the range 10~^ — 1. 

In §2 we discuss magnetized accretion flows and 
the importance of Ohmic dissipation in addition to 
the earlier considered viscous dissipation. We argue 
that the Ohmic heating is due to plasma instabilities 
which heat the electrons. In §3 we treat a model for 
the radial variation of electron and ion temperatures 
assuming that a fraction g of the dissipated energy 
goes into heating the electrons and a fraction (1 — 
g) goes into heating the ions. The electrons cool by 
bremsstrahlung and cyclotron radiation and exchange 
energy with ions by Coulomb collisions. In §4 we 
discuss conclusions of this work 



In quasi-spherical accretion onto a compact object 
of mass AJ of Schwarzschild radius = 2GAd/c? (for 
a black hole) the accreting matter is likely to be per- 
meated by a magnetic field B(r,t). Typically the 
accreting matter is ionized and consequently highly 
conducting with the result that the magnetic field 
is frozen into the flow. One result of this is that 
\Br\ oc r^^ . Thus the magnetic energy-density varies 
as £mag — B^/87r (X r~^. On the other hand the 
kinetic energy-density varies as Skin = oc r~2. 

Thus one can expect that equipartition between mag- 
netic and kinetic energy-densities occurs in the flow 
at a large distance r = r^qui ^ fg (Shvartsman 1971) 
and that it is maintained for smaller r. Further accre- 
tion for r < Tequi is possible only if magnetic flux is 
destroyed by reconnection and the magnetic energy 
£mag is disslpatcd. The dissipation of magnetic en- 
ergy was first taken into account by Bisnovatyi-Kogan 
and Ruzmaikin (1974) who showed that accretion for 
conditions of equipartition {Smag £kin) is accom- 
panied by the dissipation of magnetic energy into 
heat with entropy s (per unit mass) production rate 
pT{ds/dr) = -3BV(167rr). We point out that the 
Ohmic dissipation of the magnetic energy is an impor- 
tant, possibly dominant heating process in advection- 
dominated accretion flows with Emag ^ £kin- In this 
regard note that although Narayan and Yi (1995) as- 
sume an equipartition magnetic field, they do not con- 
sider the Ohmic heating. 

The basic equations for accretion flows with £mag ~ 
£kin are 

( 1^ + V • V ) V = --Vp + g+ —3 xB + v^V^w 
\dt J p pc 

(la) 

— = Vx (vxB)+77„v2B , (16) 

where v(r, t) is the flow velocity, p{r, t) the pressure, 
g = — V(GM/|r|) the gravitational acceleration, 
the microscopic kinematic viscosity coefficient, and 
r]m the microscopic magnetic diffusivity. 

It is well known that the microscopic classical 
transport coefficients Vm and 7]^ are much too small 
to directly influence on the macroscopic flow v and 
magnetic field B evolution. For example, for condi- 
tions pertinent to a flow onto a massive black hole, 
n ~ 10i2cm-3, Ti ~ lO^^K, and B ~ lO^'G for 
r ~ rs, the Reynolds number for the flow Rcy = 
'>'\y\/i'm ~ 10^^, where Um ~ I'li/Tu is the viscosity 
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appropriate for a tangled magnetic field (Braginskii 
1965; Paczynski 1978), and where Vgi ^ 10-^cm is the 
ion gyro-radius and m ~ lO^s is the ion-ion Coulomb 
scattering time, and ujdTu ^ 1 with ujd ^ 10*/s the 
ion cyclotron frequency. [Under some conditions it is 
possible that Vm is larger than r^i/ru as discussed by 
Subramanian, Becker, and Kafatos (1996)]. The mag- 
netic Reynolds number Rcb = r\v\/r]m ~ 10^^, where 
'Hm = c^/(47rcrs) with as the Spitzer conductivity. 

It was proposed by Shakura (1973) that accretion 
flows are in general turbulent and that roughly equa- 
tions (1) should be taken with turbulent transport 
coefficients Vt and r/t replacing the microscopic coef- 
ficients, and with v — > v and B — > B interpreted as 
mean fields. The turbulent viscosity has a crucial role 
in thin Keplerian disks where it provides a mechanism 
for the outward transport of angular momentum. 
According to Shakura (1973), vt = acgiH, where 
a = const, is the above-mentioned dimensionless vis- 
cosity parameter, Csi is the ion sound speed, and H 
is the half-thickness of the disk which is the outer- 
scale of the turbulence. Note that for an advection- 
dominated accretion flow, H ^ r. The shear stress 
in a magnetized accretion flow, which causes outflow 
of the angular momentum, appears in large part be 
due to magnetic stress (Eardley and Lightman 1975; 
Brandenberg et al. 1995; Hawley, Gammie, and Bal- 
bus 1995). Bisnovatyi-Kogan and Ruzmaikin (1976) 
argued that iji ^ vt- The turbulent diffusivity will 
have a crucial role in dissipating the magnetic energy 
in advection-dominated flows. In addition to vt and 
rjt, there will be a turbulent transport coefficient ah 
(with units of cm/s) associated with the helicity of 
the turbulence in a rotating accretion flow (see, for 
example, Ruzmaikin, Shukurov, and Sokoloff 1988). 

Neglecting for the moment the possible difference 
between Tg and Tj and the radiative energy losses, 
energy conservation for the accretion flow can be ex- 
pressed in terms of the mean fields as 

rr^ds I f_ _ 2^ _ V 

+i-%(VxB)^ (2) 

where s is the entropy per unit mass. The first term 
on the right hand side of (2) represents the viscous 
dissipation or heating of the plasma, and the sec- 
ond term the Ohmic dissipation. The two terms are 
of comparable magnitude for an accretion flow with 
Smag ^ ^kin and Vt r]t. 



However, equation (2) says nothing about the ac- 
tual microscopic dissipation of energy in the plasma. 
Rather, it expresses the loss of energy from the outer- 
scale (^ r or i7 if i7 <C r) of the flow v and from the 
B fleld by the nonlinear processes implicit in equa- 
tions (1) and the presumed Kolmogorov cascade of 
this energy to smaller scale eddies and field struc- 
tures of the flow. The turbulence may be charac- 
terized by wavenumber-frequency ensemble averaged 
spectra (v^^) and (B^^), where the wavemimber 
ranges from the small value corresponding to the men- 
tioned outer scale kmin ~ 'r~^ to some much larger 
value kmax ^ krnin- Thc Conventional Kolmogorov 
description has a dissipation scale corresponding to 
kmax ~ {Re)'*kmin which correspouds to an unphys- 
ically small length scale using either Re^ or Res- 
Thus, the actual dissipation must be due plasma in- 
stabilities. 

The relevant plasma instabilities are probably cur- 
rent driven in response to the large mean electric field, 
E = —V X B/c — a?iB/c-|-rytV x B/c, which in general 
has a significant component parallel to B. It is un- 
clear to us why current driven instabilities resulting 
from S|| were not considered by Begelman and Chi- 
nch (1988). The typical electric field |E| - 10^ V/cm 
(for r ~ rs) is much larger than the Dreicer elec- 
tric fleld for electron runaway (Parail and Pogutse 
1965), Ed = 4TTe^ine/kTe)inA - lO"" V/cm for 
Tf. ~ 10^ K, where rie is the electron density. Thus 
the electrons will runaway. An electron becomes rel- 
ativistic in a distance of travel of ~ 1 cm which is 
comparable to the electron gyro radius. The drift 
speed of the electrons parallel to B will be sufficient to 
give rise to streaming instability (Parail and Pogutse 
1965). Streaming instability will occur if the electron 
drift velocity is larger than the ion thermal speed. 
In contrast with the ions, the travel distance for a 
proton to become relativistic is ~ 10^ cm. However, 
acceleration of protons parallel to the magnetic fleld 
is strongly suppressed by scattering by magnetic fluc- 
tuations (Alfven waves) with wavelengths of the or- 
der of the proton gyro radius which are generated by 
the proton streaming (Kulsrud and Pearce 1969). For 
these reasons we believe that most of the free energy 
driving the instability goes into heating the electrons. 
Howevcir, we also consider the case where a fraction 
g of the dissipated energy goes into heating the elec- 
trons and (1 — g) goes into heating the ions. We il- 
lustrate thc behavior in this case with the following 
simple model. 
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3. Model 

Wc generalize equation (2) by taking into account 
(a) that Tj and T,. may differ with energy exchange 
between ions and electrons by Coulomb collisions, (b) 
that the Ohmic plus viscous dissipation heats elec- 
trons and ions as discussed below, and (c) that the 
main energy loss is from optically thin bremsstrahlung 
and optically thick cyclotron emission. Note that the 
thickness of the flow H/r is not restrirttxl. Note also 
that in contrast with Narayan and Yi (1995), no as- 
sumption is made that a constant fraction / of the 
dissipated energy is advected inward. Hence 

3 dTj Tj dp ,^ ,„ „ , ,„ > 



3 dTe 

2~dt' 



Tedp 
p dt 



9'H—Cbrem—Ccyc+l^ie{Ti — Te) , (36) 



where <? < 1 is the fraction of the Ohmic plus viscous 
dissipation which goes into heating the electrons. We 
assume g = const, which we view as more physically 
plausible than the assumpition that / = const, of 
Narayan and Yi. For simplicity of the formulae we 
assume Tj < mjC^ and < me(?, where Tj and Tg 
are measured in ergs. Here, 



4(27r)2ne4 / T, 



rriinie 



nie 



inA 



is the ion-electron energy exchange rate with £nA = 
0(20) the Coulomb logarithm (Spitzer 1940); H w 

(9/4)TOjQ:(c5i/?;x)^ v\J jr is the heating rate per ion 

with J7 = l-(rs/r)^; Cbr-em « naraf meC^{Te/mec'^)i 
is the bremsstrahlung cooling rate per electron with n 
the electron or ion density, <jt the Thomson cross sec- 
tion, and ttf the fine structure constant; and Ccyc ~ 
TeUi^^Ml/{8Tr^nc^r) is the self-absorbed cyclotron ra- 
diation cooling rate per electron, with Mc ^ 1 
the cut-off harmonic number of the cyclotron ra- 
diation below which the radiation is self-absorbed 
(Trubnikov 1958). For Mc > (2/9)^ > 1, with 
II = nieC^/Te, Trubnikov's analysis gives Mc ~ 
(2/i/9)(l + en{V)/iif, where V » uj%r/iaOceMc), 
with LUpe and Uce the electron plasma and cyclotron 
frequencies respectively. Trubnikov's expression for 
Ccyc is similar to that of Narayan and Yi (1995). 

It is useful to rewrite equations (3) in dimension- 
less form. Note that d/dt = Vr{d/dr) with Vj. = 

— {2>/2)aTiVKi and that H/r = , number density of 



electrons or ions n = M / {Q-namir'^'tj vk), mass den- 
sity p = nra-i , magnetic field B = \1Mvk I {"iar^Tj )] 2 , 
where Tj = Ti/Ty with Ty = GMrrii/r the virial tem- 
perature. We also normalize the electron temperature 
with the same Ty, Tg = Tg/Ty. Equations (3) become 



dr 



= -{i-g)n + A{fi-fc) , 



(4a) 



dT . , . 

= -[{2 + O9 - + Cbrem + Ccyc 

-(2 + C)i(f,-fe) , (46) 

where f = r/rs, with rs the Schwarzschild radius, 

= Tc/Tj_, and where 



7^ 2 

2^2 



H 



^brem 



J_ 

2f 




(4c) 
(4e) 



(4/) 



where Le = ATrGMrriic/ax is the Eddington lumi- 
nosity, and Te = e^/(meC^) is the classical radius 
of the electron. The terms dTi/df and dTg/df in 
equations (4) describe the advection of energy by 
the flow. Apart from the cyclotron cooling the dif- 
ferent terms depend only on a and M(?/Le- The 
cyclotron cooling is relatively more important for ac- 
cretion onto a stellar mass object than for accretion 
onto a massive black hole. The assumed condition 
for optically thin bremsstrahlung radiation requires 
{Mc^/aLE)f-^ < 1 for T- = 0(1). 

We have solved equations (4) starting from differ- 
ent given 'initial' values of T; and Tg at large f = 10^, 
different accretion rates Mc^ = (0.01 — 1)Le, values 
of a, and values of 5 = — 1, and integrating inward. 
For the accretion rates where advection-dominated 
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flows are suggested to occur (Narayan and Yi 1995), 
Mc^ < 0.1 Lb for a = 0.1, we find that the the scaled 
ion temperature Tj remains almost constant, whereas 
the scaled electron temperature Tp decreases rapidly 
as f decreases from 10^. In this limit, the Coulomb 
energy exchange between ions and electrons is negli- 
gible. The advcction terms on the left-hand-side of 
equation (3b) are also negligible. Consequently, the 
Ohmic heating of the electrons gli goes into radiation, 
mainly cyclotron radiation; that is, gH ~ Ccyc- The 
total radiation is the volume integral of gJin which 
gives gGMM/{2ri), where rj is the inner radius of 
the flow. Thus, the radiative efficiency is reduced by 
a factor of g from that of a thin disk with Tj = Te <C 1 
which is the volume integral of Tin. This efiiciency 
can be very small compared with unity only if g is 
very small compared with unity. 

4. Conclusions 

This work considers magnetized advection-domin- 
ated accretion flows where the magnetic field is in 
equipartition with the turbulent motions of the flow 
(Shvartsman 1971). The magnetic energy density of 
the flow must be dissipated by Ohmic heating with 
a rate comparable to that of the viscous dissipation 
(Bisnovatyi-Kogan and Ruzmaikin 1974). We argue 
that the Ohmic and viscous dissipation must occur 
as a result of plasma instabilities. Further, we argue 
that the instabilities arc likely to be current driven 
in response to the electric field (associated with the 
turbulent motion) which has a significant component 
parallel to the magnetic field. These instabilities are 
likely to heat mainly the electrons. We have anal- 
ysed a model for the radial variation of the electron 
and ion temperatures assuming that a constant frac- 
tion g of the viscous plus Ohmic heating goes into 
heating the electrons and a fraction (1 — 5) goes into 
heating the ions. In contrast with Narayan and Yi 
(1995), we do it not assume that a constant fraction 
/ of the dissipated energy is advected inward by the 
flow. The electrons cool by bremsstrahhmg and cy- 
clotron radiation and exchange energy with the ions 
by Coulomb collisions. At large accretion rates M, 
Coulomb collisions act to give Ti « Te, high radiative 
efiiciency, and geometrically-thin, optically-thick disk 
accretion. For small accretion rates, where advection- 
dominated accretion flows are suggeested to occur, 
and only Coulomb energy exchange between ions and 
electrons, a regime of optically thin accretion flows 
with a large difference between ion and electron tem- 



peratures {Te <C Tj) exists (Shapiro, Lightman, & 
Eardley 1976). Here, we emphasize that the accretion 
flow properties depend critically on the Ohmic heat- 
ing of the electrons. For small accretion rates where 
the electron temperature is much less than the ion 
temperature, we show that the Ohmic heating of the 
electrons gives a radiative efficiency which is reduced 
by a factor of g from that for a thin disk. Thus, the 
tiny radiative efficiencies (< 10~^) found by Narayan 
and Yi (1995) correspond to tiny values of g which 
are unlikely for the reasons discussed in §2. 

Plasma instabilities due to electron-ion streaming 
(for electron drift velocity larger than the ion ther- 
mal speed) may greatly enhance the energy exchange 
between ions and electrons. In this case the two- 
temperature regime disappears, the ion and electron 
temperatures collapse to small values, Tj^e ^ Ij and 
the disk is geometrically thin. That is, advection- 
dominatcd accretion flows do not occur (Fabian & 
Rees 1995). 
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